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ABSTRACT 

The fast, dense winds which characterize Wolf-Rayet (WR) stars obscure their underlying cores, 
and comphcate the verification of evolving core and nucleosynthesis models. Core evolution can be 
probed by measuring abundances of wind-borne nuclear processed elements, partially overcoming 
this limitation. Using ground-based mid-infrared spectroscopy and the 12.81/xr7i [Neil] emission line 
measured in four Galactic WR stars, we estimate neon abundances and compare to long-standing 
predictions from evolved-core models. For the WC star WR 121, this abundance is found to be > 11 x 
the cosmic value, in good agreement with predictions. For the three less-evolved WN stars, little neon 
enhancement above cosmic values is measured, as expected. We discuss the impact of clumping in WR 
winds on this measurement, and the promise of using metal abundance ratios to eliminate sensitivity 
to wind density and ionization structure. 

Subject headings: infrared radiation — stars: Wolf-Rayet — techniques: spectroscopic 



1. INTRODUCTION 

Wolf-Rayet (WR) stars are evolved massive stars 
characterized by high mass loss rates (IQ-'^ < M < 
10^4 M© yy-'^^MwR ^ 10^° M©), driven in fast {v^ - 
few X lO'^km/s) stellar winds. The dense WR winds ob- 
scure their underlying cores and the region from which 
the outflowing wind material is initially accelerated. The 
difficulty in quantifying bulk parameters of WR stars 
which arises due to this obscuration is typified by the 
uncertainty regarding the appropriate photospheric tem- 
perature to assign them. When the most commonly 
used definition of the photosphere is adopted (the lo- 
cation of optical depth t ~ 2/3), a temperature degen- 
eracy arises in the models of WR atmospheres — for 
all WR subtypes, the effective photospheric tempera- 
ture derived is 12/3 ~ 30, 000 K, for reasons unknown 
l)Schmutz. Lcither er. fc Grue nwald 1992). More effective 
for predicting luminosity and emergent flux distribution 
is the inferred core temperature, but this assignment de- 
pends critically on the ad hoc choice of a velocity struc- 
ture of the wind, the commonly assumed form of which 
iSchmu tz ( 1997) showed to be largely invalid in the sin- 
gle WR star (HD 50896 - WR6) for which the optically 
thin, supersonic portion of the wind has been modeled 
hydrodynamically. Testing models of WR core evolution, 
and the advanced nuclear reactions which occur there, is 
complicated by this disconnect between the observable 
bulk properties of the wind, and the wind-driving core 
buried beneath it. 

A powerful technique for probing WR core evolution 
which sidesteps these difficulties is available in the mea- 
sured abundances of wind-borne nuclear processed ele- 
ments. Neon in particular undergoes a remarkable abun- 
dance change during the later stages of a WR star's life- 
time. By the end of the more evolved WC phase (char- 
Electronic address: [j dsmith@ as. arizona. edil] 



acterized by wind material dominated by a-burning by- 
products, carbon in particular), ^^Ne becomes the fourth 
most abundant element, after ''He, "'^C, and ^*^0. The 
reactions of interest contributing to the creation of neon 
during He-burning in massive stars are (^aedcr 1983): 
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Essentially all of the '^"'N produced via the CN branch 
of the CNO cycle which dominates the earlier WN evo- 
lutionary phase is converted to ^^Ne. The further con- 
version of neon to ^'''Mg and ^^Mg is inefficient except 
at the highest temperatures (achieved only in stars with 
initial masses > lOOM©), and the production of ^°Ne via 
^^O is also negligible except in the most advanced WC 
and WO stages of the highest mass stars. The two main 
consequences of these critical neon production chains are 
a strong increase in the overall abundance of neon by a 
factor of ^200 over the course of the WR phase, and 
a rise in the isotopic abundance ratio ^^Ne/^^Ne from 
'--^O.l to ~35. Both of these changes to the neon abun- 
dance occur quite rapidly (in the course of several thou- 
sand years) at the onset of the WC phase. Late-stage 
depletion of ^^Ne by conversion to magnesium is only 
30% during the final 10^ yr of the most massive stars' 
lives, with the overwhelming majority of WRs maintain- 
ing their full neon excess, thanks to the combined effects 
of the interior mixing and mass-loss which bring material 
to the surface. Assuming normal p-p nuclear processing 
of hydrogen entirely to helium, the cosmic abundance of 
neon by number is Ne/He = 3.7 x lO"* (see S l^ . The 
neon abundance predicted by WR evolutionary models 
is Ne/He — 6.6 x 10"'^, or over 17 ti mes the cosmic 
value. This result, first described by [Maeded fl983), 
has remained valid despite recent model updates to acco- 
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modate rotational mixing and the tu rbulent diffusion of 
core material into the wind it drives l)Maeder fc Mevneil 

We will give background on the detection of neon in 
WR winds in §|21 describe the observations in §|21 out- 
line the abundance calculation and inputs in §§0H3 and 
present the results in §|S1 

2. BACKGROUND 

Despite lingering uncertainties concerning the true 
structure of WR winds, broad agreement between the 
surface abundances predicted by core evolution models 
was obtained earl y on for almost all a bundant elements, 
both for WN fe.g. lCrowther. Hillier "fc Smi th 1995a) and 
WC (e.g. Macder fc Mevneti 119941) stars. However, a 
long standing discrepancy concerning the model-sensitive 
neon abundance predictions centered on 7^ Velorum 
(WC8), the nearby, optically-brightest WR star. While 
the Maeder abundances were well-confirmed for all other 
elements, the measured neon abundances of 7^ Vel re- 
mained quite low, near or just above the cosmic value. 

In one of the ear liest mid-infrared (MIR ) mea sure- 
ments of a WR star, lAitken. Roche, fc Allenl (|1982D cal- 
culated Ne"*" an d S'^"'" abundances in 7^ Vel from ground- 
based spectra. Ivan der Hucht fc OlnoiJ l)1985j) found a 
similarly elevated neon abundance for 7^ Vel using IRAS 
LWS spectra, seemingly confirming its predicted over- 
abundance with respect to cosmic levels. Both of these 
results, however, w ere later significantly corrected by 
iBarlow et all l)1988() . who discovered flaws in the cal- 
culations and in the atomic inputs, and found a re- 
vised neon ab u ndanc e again quite close to cosmic values. 
iDessart et al.l l)2000 l) adjusted 7^ Velorum's neon abun- 
dance yet again, using ISO data and an improved dis- 
tance value along with more modern, dumping-corrected 
mass loss rates to derive a value c oincidentally quite 
close to the original determination of lAitken et al.L and 
in good agreement with theory. IWihi^t alJ l|199^ 
found an elevated neon abundance in the ISO SWS spec- 
trum oOZ2JiiiWC5]i_andJnJotal, four ISO WC 
stars llWillis et al.ll997tlDessart et al.i.2nom and one WN 
star llMorris et al j200(H) have vielded neon abundances. 
iMorris e't'ldTl(2004tl used early Spitzer IRS spectra to 
measure the neon abundance of WN4 star WR6, and 
found values consistent with no enhancement. 

3. OBSERVATIONS 

'S mith fc Hoild^ lf2nnl (hereafter!SHQl') present a flux- 
calibrated 8-13 /xm spectral survey of a large sample of 
northern Galactic WR stars representing all sub- types. 
Among the sample, four stars exhibited broad, non- 
nebular [Neil] 12.81 /im line emission. These are listed 
in Table ^ along with spectral types, photometric or 
cluster distances and reddening, visual magnitudes and 
mid-infr ared fluxes. The spectra were obtaine d with 
Score l|Smith et a,].ll199SllVa,n Gleve et a,lJIT99a. a pro- 
totype instrument for the short wavelength, high resolu- 
tion module of Spitzer 's IRS spectrograph f Houck et al.l 
Hq^), operated at the Palomar 5m^. Standard beam- 
switched SHz chopping and nodding were used to re- 

^ Observations at the Palomar Observatory were made as part 
of a collaborative agreement between the California Institute of 
Technology, the Jet Propulsion Laboratory and Cornell University. 
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Fig. 1. — Score spectra of two of the program stars. The [Nell] 
12.81 fim line is inset for all four stars (solid: WR105, dotted: 
WR116, see SHOl for full spectra). The dusty WC9d has a steep 
slope typical of ~1000K heated dust emission, with a slight break 
in the slope at 9.5 fj,m, characteristic of moderate, broad silicate 
absorption. The noise excess in both spectra at 9.7/im is due to 
telluric O3. 



move the sky signal. The chop amplitude was cho- 
sen to be small enough so that the object fell within 
Score's 12" diameter slit-viewer field when the slit was 
on adjacent sky. Two equal amplitude source images 
in the slit-viewer's 11.3/im silicate filter were thus ob- 
tained simultaneously with the spectra, and used to cor- 
rect the absolute flux calibration for the changing slit 
throughput function, which is affected by seeing, point- 
ing accuracy, and object acquisition. This calibration, 
and a general description of the spectral reduction pro- 
cess, are described in more detail by SHOl. Initial flux- 
calibration was performed using same-airmass observa- 
tions of Gohen et al. (1995) infrared standard stars. Line 
strengths were computed using a polynomial fit to the 
nearby continuum. Fig.^shows the Score MIR spectra 
of one of the WN and the single WG program stars. The 
WN star exhibits lines of helium and neon, while the WG 
star exhibits only [Ne 11] . 

One distinct difference between the Score data used 
here and ISO SWS spectra used to compute previous 
neon abundances is worth mentioning. The entirety of 
ISO's SWS spectral beam (~ 20" x 35") is mapped 
onto the spectrograph's large detector elements and is 
included in the recorded spectra. For this reason, con- 
tamination by nebular emission lines originating in the 
often bright, extended nebulae surrounding WR stars is 
often a concern. Score's small 1" x 2" slit excludes most 
of this nebular emission, leading in some cases to signifi- 
cant differences between spectra of the same non- variable 
W R star observed wit h bo th ISO and Score (cf. WR146 
in lWillis et all l[T99l vs. ISHfill as described in §4.2 of 
the latter). While we know of no cases in which neon 
or other abundances computed using ISO data were af- 
fected by this type of contamination, the SCORE spectra 
of fainter WR stars used here should be relatively less 
affected by nebular emission. 
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TABLE 1 

Program WR Stars with Observed Neon 



Object 


Name 


Spectral 


D 


Ay 


V 








Type 


(kpc) 


(mag) 


(mag) 


(Jy) 


WR 105 


Ve2-47 


WN9h 


1.58 


8.63 ±.03 


12.92 


0.68 


WR116 


STl 


WN8h 


2.48 


6.89 


13.38 


0.44 


WR121 


AS 320 


WC9d 


1.83 


5.72 ±.32 


12.41 


1.51 


WR124 


209 BAG 


WN8h 


3.36 


4.43 


11.58 


0.14 



N ote. — All spectral types, magnitudes, and extinction coefficients are 
from Ivan de r Hucht (2001). Extinctions with quoted errors are for stars 
with known cluster associations, and are considerably more accurate. 




Fig. 2. — Wind model geometry. Observer is to the right, along 
the line of sight coordinate x. The clumps of fill factor 5 are indi- 
cated schematically. The shading to the left of the core along the 
line of sight indicates the region from which no continuum emis- 
sion is observed. Also indicated are representative physical radii for 
the c ore, sonic point, and the photosphere (from Heeer &: Lanee'3 
I1996r , as well as the location of the continuum optical depth = 1 
at a v ariety of wavelengths, from the WN5 wind model of .Hillier. 
USHl. 



r(6) = J SKffdx = SIC{iy,T) J n^ndx — ^ 26h^ ^^'^ 

(3) 

where the reduced free-free opacity IC{iy,T) = Kff/nue, 
7e is the number of electrons per ion, and we have 
made use of the l/r^ density profile of Eq. |21 As- 
suming a constant, thermal source function, and inte- 
grating over cylinders of constant impact parameter b 
(and thus constan t free- free optical depth), we recover 
I Wright fc Barlo'wi s infrared/radio free- free flux expres- 
sion, modified to include the effects of clumping via the 
fill factor 6: 



MZ 



4/3 , -2/3 
I £)2 



2/3 



c mym. 



3/2 



(4) 



where D is the distance to the star, gi, is the frequency de- 
pendent free-free Gaunt factor, and Z is the rms average 
charge per ion. Most WR mass loss rate estimates are 
derived from radio measurements of the free-free emis- 
sion using Eq. 0| Given the same assumptions of atomic 
parameters of the wind [Z , /z, 7e), it is apparent that, in 
the absence of information about the clumping fill factor 
5, the clumping -scaled mass loss rate, Msd = M/VS, is 
derived. 



4. CALCULATING NEON ABUNDANCE 

4.1. Mass Loss & Clumping 

For a uniform, spherical, but clumped wind with ter- 
minal velocity Voo and constant volume filling fraction S 
{0 < 6 < 1), as depicted in Fig. [21 the mass loss rate can 
be written: 



M — Airr 6 ■njimHVo 



(1) 



where /i is the mean atomic mass per ion, and n is the 
number density of the ionized gas. Defining the standard 

mass loss coefficient A = -, — — , the density can be 

expressed as: 



Jj,2 



(2) 



The dominant radiative output of the wind at mid- 
infrared and longer wa velengths is free-free emission 
l|Wright fc Barlow n'975ft . The free- free optical depth 
along a particular line of sight through the clumped wind 
to the observer is 



4.2. Two level emission 

For fine structure lines arising from ions with ground 
states consisting of only two energy levels, the frac- 
tional abundance of the ion by number, ji = Uij rij, 
can be calculated straightforwardly from the observed 
line flux by neglecting all other transitions. Following 
iBarlow eTa l. (1988), the flux due to a given line transi- 
tion can be written: 



(5) 



where A^i is the Einstein emission coefficient for the line 
in question, and n„ is the density of the ions populating 
the upper level of the transition. 

The upper level density and electron density follow 
from the definition of 7^: 



]uiiA 



leA 



(6) 



4 



Neon Abundance of Galactic WR Stars 



where /„ is the fraction of that ionic species present in 
the upper level. Plugging into Eq. I^for wc find the 
line flux to be: 

oo 

Tui = jy^ / fu{ne{r),T) dr (7) 



To derive the upper-level fraction fy^{n^{r),T), we be- 
gin with t he total coUisiona l de-excitation rate per unit 
volume is (|Osterbrocklll974) : 



nenuQui 



2tt flui{T)nenu 




1712 



(8) 



where ftui {T) is the collision strength, and tj„ is the sta- 
tistical weight of the upper level. The coUisional exci- 
tation rate can be obtained from the de-excitation rate 
using detailed balance: 



neTiiqiu 




2tt fliu{T)neni 



(9) 



In the simple two level ion, statistical equilibrium be- 
tween the levels can be written succinctly: 



neUiqiu = neTiuqui + 



(10) 



from which the upper level population can be found sim- 
ply, using m = Uu +ni, as 



1 + 

Qui , 



(11) 



where ric = is the critical density in the two level 
approximation. Plugging in for rig from Eq. |S1 and the 
ratio of the coUisional rates from Eqs.|Hl-|ni we find 



£)2 



1 + 



UcSr^ 

leA 



-ekT 



dr 



(12) 



from which we can solve for the fractional abundance by 
number of the ion contributing to the line emission: 




(13) 



Notice that 7,; (x ^fE/A^ . In the absence of informa- 
tion about the clumping fill factor 5, the scaled mass 
loss rate is derived from radio continuum measurements 
(see S I4.1II . If the scaled mass loss rate is known then 
A c<. M = MsciVS, and the weak scaling of fractional 
abundance with clumping fac tor becomes 7i oc 

Note that we di ffer from iDessart et al.l ()200C1|) and 
iMorris et al.l 1)200(1 ). who perform a numerical integra- 
tion over the upper level ionic fraction /„, after trans- 
forming from radial to density coordinates to mitigate 
finite step size inaccuracies near the origin. While they 
find discrepancies of order 20% between their analytical 
integral analogous to Eq. 1121 and their numerical inte- 
gration, we find no reason that such differences should 
occur for this two-level transition (for which the upper 
level fraction integral is exact) and presume it must have 
arisen from different atomic data inputs in the statistical 



equilibrium code used, or the finite numerical resolution 
of the integration. 

Typically, WR wind abundances are formulated with 
respect to helium, almost always the most abundant el- 
ement in all but the least evolved, late-type WN winds. 
Given knowledge of other elemental abundances, this can 
be calculated according to: 



Xi_ _ / iL ii 

Hc^^'r^Ifc^Hc^He'^He 



(14) 



where represents the ion i, and the standard short- 
hand X/Y = Nx/Ny has been used. Only the most im- 
portant additional elements beyond helium are shown, 
and in many cases, the abundances of most will be so 
low that only one other element besides helium need be 
considered (e.g. carbon for WC stars, nitrogen and/or 
hydrogen for WN stars). 

5. INPUTS 

The fractional ionic abundance can be calculated di- 
rectly from Eq. ^1 given estimates of the distance, 
atomic wind parameters, and reddening. Table |21 lists 
the atomic data used for the two fine-structure lines con- 
sidered here. The measured neon and sulfur line fiuxes 
used here differ slightly from those presented in ..SHOll . 
though the reduced spectra are unchanged. This is a 
result of better estimation of the continuum. 

5.1. Chemical Composition 

The adopted chemical composition of the winds serves 
only to normalize the computed abundance with respect 
to helium (Eq. 114(1 , but does not change the results oth- 
erwise. All abundances mentioned are by number. 

The three WN stars with neon present are all late 
types — the only WR type with any significant hy- 
drogen remaining. The abund ances assumed follow 
NuEis. Crowth er. fc Wil lis (1998), with N/He=0.005 for 
aU late WN's. The hydrogen content (H/He) of 
WR 105 (2.3) a n d WR 124 (1.9) were available from 
INugis fc LamersI l|2000D . computed using optical Hei, 
Hen and Hi line measurements. The latter value is also 
adopted for the WN8h star WR 116. 

For the carbo n star WR121, the weighted mean in 
iNugis fc LamersI for subtypes WC8-9 of carbon abun- 
dance (C/He=0.18), as well as oxygen abundance 
(O/C=0.2) was used. The WC9's hydrogen abundance 
was assumed to be zero. 

5.2. Mass Loss and Terminal Velocity 

The mass loss rate enters Ea.ll3lthrough the mass loss 
coefficient as yl'^/^, and hence significantly affects the 
values obtained for the abundance. Unfortunately, rates 
derived for individual stars often differ substantially, de- 
pending on the input assumptions, and details of the 
measurement. Along with imprecise distance estimates, 
poorly-constrained mass loss rates introduce the largest 
uncertainties in the computed abundances. 

Eq. ^ can be used to determine the mass loss rate 
(or at least the scaled rate — Mj^fS) from the mea- 
sured radio flux, given independent estimates of the dis- 
tance, terminal expan sion velocity, and atomic pa rame- 
ters of the wind (e.g. iLeitherer et all 119951 11997)1 . The 
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TABLE 2 

Atomic Data for Neon and Sulfur Lines. 



Ion Transition Wavelength oj; uju f^ul ^ul ric (8000K) 

(fim) (lO^^s'i) (cm-^) 

[Siv] ^^"3/2 ^ ^^1/2 10.5105 2 4 8.47 7.70 3.77x10* 

[Neil] ^J'1/2^^^3/2 12.8136 4 2 0.28 8.59 6.36x10^ 

Note. — n^i computed at Te = 8000K. IDessart etld] ^KKt) . 



radio emitting regime is quite far out in the wind, such 
that the assumption that the terminal velocity has been 
obtained is always valid. The infrared free-free contin- 
uum, however, may not be optically thick (recall that 
K// oc in the long wavelength limit). Contributions to 
the 10 /xm continuum from the underlying stellar photo- 
spher e are typically only 10% of the free-free wind contin- 
uum l|Barlow et al.ll98l|) . and the cross-over from photo- 
sphere to wind-dominated emission occurs near 1-5 /xm; 
however, since the free-free opacity is lower, the infrared 
emission originates deeper in the wind, sampling out- 
flow which may not yet have reached terminal velocity, 
such that the assumption leading to Eq. El and a 1/r^ 
density distribution is violated, and the spectral index 
steepens compared to the constant velocity limit. ^ In 
the absence of a valid velocity law, one possible tech- 
nique to measure the mass loss rate is to scale the in- 
frared flux with an empirically determined spectral index 
a (fu oc zy"), as iBarlow et alJ l|1981|) did for 7 Velorum 
and HD 192163 to estimate a scaling with a = 0.76. 
More recently, Crowthcr et al. ( 1995a) found a value in 
close agreement, a = 0.74, based on a tailored analysis 
of several WN stars. 

Scaling our MIR fluxes with this method and spectral 
index a — 0.74, mass loss rates for the three program WN 
stars (i.e. those without clear evidence of dust emission) 
were calculated. Each spectrum was dereddened using 
a composite extinction curve of Cohen (priv. communi- 
cation), which, longward of 4.7 jum, joins smoothly the 
curves used bv iCohenI l|1993D . The small inferred MIR 
reddening corrections are relatively uncertain, but the 
largest uncertainty is for the dusty types not considered 
in this calculation. 

A 0.2 /xm wide region centered on 12 /xm was averaged 
to arrive at the SCORE flux density. IRAS Point Source 
Catalog 12 /im flux densities (Cohen 1995) were found to 
be 2-5 X higher, likely due to the very large IRAS beam- 
size. The extrapolation of the de-reddened 12 /im fluxes 
was performed to 4. 8 GHz (6.25 cm) for comparability 
with iLeitherer et al.l il997D. Atomic pa rameters of the 
wind were adopted from ILeitherer et aO . using the usual 
relation for mean ionic mass 

A*=^F^ (15) 

where rrii is the atomic mass of the ion. Only H, He, 
and C are included, since O never reaches abundances 

^ A simple understanding of why the spectral index must steepen 
in an accelerating wind is had by noting that at different frequen- 
cies, T"// ~ 1 occurs at different depths, since ^ = Kff is quadratic 
in the density. If the rate of change of density with radius is steeper 
(as in a region of acceleration), adjacent frequencies will originate 
further apart in the flow, exaggerating the difference. 



TABLE 3 

Example Mass-Loss Rates from Various Methods. 



Method log{Af) (Mo/yr) Reference 



WR 105 (WN9h) 



Radio 


<-4.41 


(1) 


Radio+Clumping 


-4.55 


(2) 


Extrapolated IR 


-4.95 


(4) 


Optical+1 /xm Line Analysis 


-4.1 


(5) 


Optical+UV Line Analysis 


-4.2 


(6) 


WR124 


(WN8h) 




Optical Line+Clumping 


-4.61 


(2) 


Extrapolated IR 


-4.7 


(3) 


Extrapolated IR 


-4.95 


(4) 


Optical+UV Line Analysis 


-3.8 


(6) 



References. — (1) ILeitherer et al.l 119971) . (2) 

INugis fc Lamer^ |200C|). (31 IB arlow et al.l 119811). (4) this 

work. (5) [achmutz^ Ham anm fc Wcss(3^wskir ~ Tl98gi) . (6) 
IHamann fc Koesterkd 119981) 



which would influen ce ii significantl y. The mean molec- 
ular weights per ion'Leithe rer et al.l compute are consis- 
tent enough within subtypes that we adopt /i — 2.0 for 
WN types later than WN6. 

The mean number of electrons per ion, 7e, and the 
RMS ionic charge, Z, were computed assuming He"*" 
is the most prevalent ion in the radio emitting region. 
These quantities are very insensitive to any other as- 
sumption, since singly ionized helium and hydrogen con- 
tribute to both cquivalently. Carbon is often assumed 
to exist as C+^ in the radio regime, but it does not sig- 
nificantly impact the calculation of either 7e or each 
of which typically only ranges from 1.0-1.2. Given the 
magnitude of other uncertainties in the calculations, we 
adopt Z = 7e = 1.1. The Gaunt factor, which enters 
Eq. ^ through the free- free opacity, depends logarithmi- 
cally on temperature and ionic charge at 4.90 GHz . We 
again follow Lcithcrcr et al. in adopting (74 s = 5.0, which 
they computed for = 10, 000 K. 

For two of the sources, mass loss rates were available 
from the dumping-corrected emission line-fitting results 
of INugis fc Lamer j (1200? ?): log(M) = -4.55 (WR105) 
and log(M) = -4.61 (WR124). Illustrating the uncer- 
tainty in the rates, TableEllists different estimates based 
on different techniques for two of the WN program stars. 
Rates based on UV, optical, IR, and radio data, with and 
without clumping accounted for, are listed for two of our 
program stars. It should be pointed out that only the 
Radio and Extrapolated IR methods do not require de- 
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tailed modeling with the implicit ad hoc assumption for 
the form of the velocity field, although the IR rates do 
depend sensitively on the measured spectral slope for ex- 
trapolation. The variance among the different estimates 
is quite large, and, in the case of WR124, the smallest 
and largest rates found differ by a factor of ~14. 

WR 105 is a confirmed non-thermal emitter 

ijChapman et alJflQQ gl). with a spectral index a = —0.3 
over the 3-6 cm radio band. Since the infrared fiux will 
likely not be modified by the thermal X-ray and non- 
thermal radio emission arising in the binary wind-wind 
collision al shocks thought to u nderly non-thermal WR 
sources (|Eichler fc Usovl I1993D . we expect this star's 
12 /im extrapolated mass loss rate (log(Af) = —4.95) to 
be more accurate than radio-based results, and adopt it 
for abundance calculations. 

Where available, we have preferred the radio and 
extrapolated infrared rates over others. In one case 
(WR116) no radio or clumping corrected results were 
available, and we therefore adopted log(M) — —4.6, 
b ased on observations o f other WN8 stars. 

iLeitherer et al.l l)1997f) put an upper limit on the mass 
loss rate of WR121 from Australia Telescope Com- 
pact Array radio no n-detection of log(M) < —4.5, and 
iBieging et al.l l)1982D fin d a similar lim i t (<- 4.55) using 
the Very Large Array. lAbbott et alJ l|1986D . however, 
had previously placed a firmer limit of log(M) < —4.8. 
We adopt log(Af) = —4.9, based on values derived from 
radio detections of other WC9 stars, and in line with the 
lower mass-loss rate WC9' s exhibit compared to earlier 
WCs ijLeithere r ct al. 1997*) . The rate inferred from scal- 
ing WR121's MIR flux is not valid, having been unduly 
influenced by the excess heated dust emission evident in 
WR121's spectrum. 

All terminal velocities are from Ivan der HuchHll200l . 
except for WR105, for w hich we adopt the value of 
lEenens fc Williams! l|1994D . Voo = 1200 km/s, which 
closely matches our measurements of the [S iv] 10.5 /zm 
FWHM in our spectrum, but deviates from most the re- 
cently cataloged value of 700 km/s. 

5.3. Temperature and Atomic Parameters 

Despite the very high inferred effective temperatures 
(>60kK) of the underlying star, WR winds are quite 
efficiently cooled by line radi ation t o < lOkK in the 
radio emitting regime (Hillier 1989). Though the exact 
temperature in the line-emitting region is uncertain, the 
temperature dependence of the derived ionic abundance 
is extremely weak at these long wavelengths {hv/kTf. ^ 
.1 =^ 7, oc T°-°^); we therefore assume = 8000K for 
all four sources. 

The mean ionic masses (//) and mean number ofelec- 
trons per ion (7e) were taken from Leitherer e t al.l l|1997f) 
in the case of WR 105 and WR 121. For WR 124, was 
computed dire ctly from the abundanc e ratios Nh/Nho, 
and Nz/Nhc of lNugis fc LamersI pOOO^ . This same value 
was used for the other WN8h star considered, WR116, 
and in both cases 7e = 1.1 was assumed from analogy 
with other late WN's. 

5.4. Total Neon Abundance 

Estimating the total neon abundance from the spectra 
is difficult, especially since the [Neiii] line at 15.56 /im 



is outside of the N-Band atmospheric window, and thus 
cannot be used to probe the gas in a higher ionization 
state. Theoretical neon ionization structure predictions 
could in principle be used to infer total neon abundances. 
Neon is now routinely included in the blanketed, non- 
LTE WR atmosphere code CMFGEN ( HilUe rfc Milled 
Il998j) . and hence its ionization structure can be modeled. 
It is an impurity species, and generally only has a small 
effect on the spectral energy distribution; it can, however, 
help drive the wind, especially in the outermost regions. 
Although a wind model specific to the very late WC9 star 
WR 121 was not available, by analogy to other models for 
low ionization WR stars (e.g. WNIO), Ne"*" is likely to 
be the dominant ion in the outer wind (J. Hillier, 2004, 
private communication) . 

Some handle on the most likely ionization state can 
also be had by noting that the ionization potential of 8"^+ 
is 35 eV, vs. Ne+: 21.6 eV and Ne++: 41.0 eV. Since 
their critical densities are reasonably close at these tem- 
peratures (ne ~ 10^), it is expected that the detection of 
[Siv] at 10.5 /xm is an excellent pre dictor of [Nei ii]. In- 
deed, in the hot WN8 star WR 147. iMorris etal] |2000) 
found [Neiii] and [Siv], but no [Neii], and non-LTE 
model predictions of the outer winds of la te WC stars 
show Ne"*" as the dominant ion species llWiUis et alJ 
IT997t) . While the presence of [S iv] clearly does not imply 
a dearth of Ne"*" (cf. WR 105, SHOl), its absence places a 
stronger limit on the possible existence of Ne"'"'' .'^ By this 
argument Ne"'" accounts for nearly all the neon in those 
stars with no [S iv] detection. Additional support for this 
conculsion as it pertains to WR121 is provided by an- 
other WC9 star recently observed with the Spitzer's IRS 
spectrograph, which offers similar resolution as Score 
but covers 10-40 /J,m. The IRS spectrum shows strong 
[Neil] without any detectable 15.56 /im [Neiii]. This and 
other early IRS WR results will be presented in a forth- 
coming paper. 

5.5. Cosmic Neon Abundance 

The "cosmic" abundance of neon measured from so- 
lar coronal lines has been revised many times over the 
past 30 years. Recommended values for the fractional 
abundance by number, ejve — 10^^ x Nnc/Nh have 
ranged from log(eNo) = 7.64 to 8.04 (c onverted from 
the total mass ab undance of, e.g ., fCamer on (1973) and 
the references in 'Maeder (1983^) . More recent values 
include those suggested by Grevesse & Sauval (199^ 
(log(eAre) = 8.08), and the value foun d from the updated 
oxygen abundances of lAsplund et alJ (^2004) (log(eAre) = 
7.84). 

For evolved WC stars, the cosmic abundance is of lit- 
tle interest, since the bulk of the neon atoms entrained 
in the wind material were created directly from helium 
burning in the core. Final neon abundances with respect 
to helium can then be compared directly to model pre- 
dictions, with little sensitivity on the initial amount of 
neon. 

The less evolved WN stars do not produce any neon 
directly, and are therefore expected to exhibit Ne/He 

Since sulfur abundance is not enhanced by nucleosynthesis and 
neon is, this argument is weaker if the Ne/S abundance ratio were 
so large that S'^"'", while present, remained undetected, despite the 
strong measured [Ne ll] emission. 
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matching cosmic abundances. However, all WR stars 
produce, and potentially consume, helium, so that the 
"cosmic" value of Ne/He actually varies, depending on 
what assumptions are used to correct the helium abun- 
dance for its enhancement a nd/or depletion via nu- 
clear proce ssing. For example, iMorris et a l.l J200QD and 
Im orris et a l. (2004) derive a corrected Ne/He abundance 
by assuming complete conversion from AH —>■ He for the 
WN targets considered: 



Ne\ 



N 



No 



Ne/H 



NHe + Nh/4 He/H + 1/4 



(16) 



Using the latest lAsplund et all abundances, this yields a 
cosmic abundance Ne/He = 2 x 10 ~^ in the severely H- 
depleted environment of WN winds. iBarlow 
perform a similar adjustment for the WC8 star 7^ Velo- 
rum, assuming C/He=0.2 to correct Ne/He for process- 
ing of all hydrogen to helium, and subsequent conversion 
of helium to further elements. 

Since any unprocessed hydrogen in the core or wind 
material (common in late WN stars) will invalidate the 
complete conversion assumption of Eq. 1161 we adopt a 
range of cosmic Ne/He using the updated solar abun- 
dances of lAsplund et al.ll|2004() . from Ne/He = 2.0x10-"* 
(complete H processing) to 3.7 x 10"* (no H processing). 

6. RESULTS 

A summary of the inputs and results of the neon 
and sulfur abundance calculations for the four program 
sources are given in Table 0] Immediately apparent is 
the overabundance of Ne+ in WR 121, the late WC star. 
The three WN stars for which Ne^ abundance was mea- 
sured are not expected to display abundance enhance- 
ments, since the byproducts of CNO processing should 
dominate. All are reasonably close to the cosmic value 
for Ne/He, ranging from l.l~4.4x cosmic (the latter cor- 
responding to the fully- processed cosmic value), with 
WR105 noticeably higher than the others. 

Presuming all neon is accounted for in Ne^, the to- 
tal neon abundance can still be increased by introducing 
the clumping fill factor S. li 6 — 0.1, the abundances 
are increased by ~1.8, presuming all other background 
abundances are unaffected. If all abundances (including 
helium) scale the same with clumping factor, the abun- 
dance ratios should technically be insensitive to it, but 
differing emission regimes complicate this argument (see 

Another interesting constraint is available by con- 
sidering the mass loss rate dependence. If for, e.g., 
the nitrogen star WR 124, the Optical+UV line anal- 
ysis based rate of iHamann fc Koesterkd l)1998|l {M = 
IO^^-^Mq yr"^) had been adopted with all other in- 
puts unchanged, the derived abundance would drop to 
Ne+/He = 2 x 10"^, or only ~5% of the cosmic neon 
value. Since the reactions which convert neon to magne- 
sium occur so rarely, and at core temperature attained 
only at latest stages of WR evolution and for the most 
massive WR cores, this casts doubt on such a high mass 
loss measurement; once predicted wind abundances are 
confirmed, neon and other products of nuclear processing 
can be used to constrain mass loss rates in the absence 
of other measures. 



The WC star WR121's neon abundance Ne+/He = 
4.1 X lO"'^, is in quite good agreement with the long- 
standing prediction of Ne/He = 6.6 x 10"'^. Though 
any significant abundance of Ne^"*", which we could not 
detect, would serve to increase Ne/He, doubly ionized 
neon is not expected to be abundant in this source, due to 
the total absence of [S iv] emission. A realistic clumping 
factor of S = 0.15 would imply an abundance exactly 
equal to the predicted value. 

Sulfur, though produced in the advanced oxygen burn- 
ing leading up to supernovae, is not significantly en- 
hanced by normal stellar nucleosynthesis, and exhibits 
no abundance increase. 



7. DISCUSSION 
7.1. Neon Detection Frequency 

Of the 29 WR stars cataloged by ISHOU only four 
exhibited measurable [Neil] emission. The absence of 
neon in the remaining 25 is likely due to two factors: 
higher typical ionization of the wind in the neon-emitting 
regime, and bright continuum from heated dust diluting 
any spectral lines which otherwise would be present. 

Among the ISO WR stars with detected neon emis- 
sion, only 7^ Vel (WC8) showed any non- nebular [Neii]. 
The remaining WC5-WC8 and WN8 ISO spectra exhib- 
ited only [Neiii] 15.55 /im, which is inaccessible in the 
ground-based SCORE data. That the low [Neil] detec- 
tion frequency is a result of higher average ionization is 
supported by the late subtypes of the four program stars 
discussed here; for both WN and WC, the subtype se- 
quence is primarily one of ionization, with early types 
exhibiting lines of increasingly higher ionization species. 
This same earlier-type, increasing ionization sequence is 
also found in the tren ds of he lium line emission strengths 
in the MIR spectra of SHGl. 

The majority of WC9 stars are dust emitters, with in- 
frared excesses 10 x or more the normal free-free excess 
seen in non-dust producing WR stars ij Williams et alJ 
Il987|) . The bright infrared continuum of heated carbon 
dust is usually unaccompanied by line emission in the 
thermal IR; presumably the lines suffer so much con- 
tinuum dilution they are unmeasu rable. Among the 6 
WC9d late-type WC's observed bv lSHOl only WR121 
showed any line emission, and indeed it is the only 
WC9 for which any MIR line emission has ever been 
observed (the two neon -emission WCIO stars reported 
by Aitke n et al.l l)1980l) based on early low-resolution 
ground-based spectroscopy were later revealed be plane- 
tary and proto-planetary nebulae). 

A notable star in the original MIR sample of 
ISmith & Houc^ without [Neii] present is WR145, a 
WN/WC type. The rare W N/WC transition stars eight 
of which are cataloged by Ivan der Huchd (|2003), 

ex- 
hibit spectral signatures intermediate between the WN 
and WC classes, and are hypothesized to represent a 
transitional stage between CNO-dominated hydrogen- 
free WN stars, and g-burning d ominated WC stars 
(jCrowther. Smith, fc Willislll995b() . A confirmation of 
enhanced neon abundance in one of the transition objects 
would support this interpretation, but unfortunately any 
neon in WR145 must be present as Ne++, as expected 
from the strong [S iv] emission the star exhibits. 
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TABLE 4 
Neon and Sulfur Abundances. 



Object 


Type 


D 


7e 






log(M) 


Line Flux 


Xi 

He , 






(kpc) 






(km/s) 


(Mo/yr) 


(10-19 W/cm2) 


(xlO-4) 




Ne+ 


WR105 


WN9h 


1.58 


1.0 


2.6 


1200 


-4.55 


1.5 


8.8 


WR116 


WN8h 


2.48 


1.1 


2.0 


800 


-4.6 


1.6 


5.3 


WR121 


WC9d 


1.83 


1.1 


4.7 


1100 


-4.9'!' 


2.1 


41 


WR124 


WN8h 


3.36 


1.1 


2.0 


710 


-4.61 


0.5 


4.2 


S3+ 


WR 105 


WN9h 


1.58 


1.0 


2.6 


1200 


-4.55 


1.2t 


0.5 



tA radio limit of log(M) < —4.8 is given by lAbbott et ahl 119861) . and we adopt the 
plausible value -4.9 based on other measured WC9 rates (e.g. WR80). 
-^Including 15% correction for helium contamination of the line. 



7.2. Neon/Sulfur Abundance 

As pointed out by iDessart et al.l 1)200(11) . the use of 

neon abundances as derived in the previous section to 
constrain true core-evolution nuclear processes is compli- 
cated by the different emitting regimes of the elements 
in question. Because of their comparatively small ra- 
diative A coefhcients, the fine structure lines originate 
further out in the wind than the helium emission, at 
densities near 10^ cm^^. A dependence on precise mass 
loss rates and (less problematic) terminal velocities also 
reduces the accuracy of direct abundance measurements, 
especially given the tendency of clumping to modify the 
measured M by factors of 2-5. A method which over- 
comes these difficulties is available in the Ne/S abun- 
dance. Since both [Siv] and [Neil] originate in roughly 
the same region of the wind, the abundance ratio derived 
from them should be independent of details of the wind 
structure, and remove any attendant systematic errors in 
estimating the mass loss rate, distance, and bulk param- 
eters of the star. 

The He i and He ii lines blended with [S iv] 10.5 /zm are 
known by wind ionization models to con tribute only very 
weakly to this line in 7^ Velorum (<15%. liDe Marco et al.l 
I2000D . Due to the blend of neutral and ionized helium 
and hydrogen lines at this location, this contamination 
fraction is expected to remain nearly constant with WR 
subtype. We have therefore reduced the measured [S iv] 
input flux in Table 0] by this amount. 

The S^+/He abundance found for the WN9 star 
WR 105 is in good agreement with the cosmic value of 
S/Hc ~ 7.5 X 10^^, and if S++ constitutes about one- 
third of the sulfur, as in WR146 ( Willis ct al. 1997|), the 
total sulfur abundance comes quite close to cosmic. The 
neon to sulfur ratio we find is Ne+ZS^"*" = 15.7, some- 
what larger than the cosmic value of the full abundance 
ratio Ne/S = 7. If, as expected, sulfur is less completely 



accounted for by S'^^ than neon is by Ne+, the total 
sulfur abundance would be increase by a larger factor 
than would neon (around a factor of 2), and the final 
implied Ne/S can approach the expected cosmic value, 
as expected in this less-evolved WN star. 

8. CONCLUSIONS 

We have computed neon abundances for 3 WN and 
I WC star from ground-based spectra of the [Neil] 
12.81 /im emission line. WR121, the WC star, shows 
elevated Ne^ abundance, with an estimated total neon 
11. Ix the cosmic value, close to the expected 17. 8x in- 
crease predicted by WR core evolution models. A sig- 
nificant population of emitting neon ions in the Ne++ 
ionization level, or realistic clumping fill factors 5 would 
each revise this value upwards by up to a factor of 2, 
but distance and mass-loss rate uncertainties also con- 
tribute. Though elevated neon abundances have been 
found in several other WC stars, this is the least-evolved 
star for which such enhancement has been demonstrated. 
The WN stars were found to have abundances close to 
cosmic, consistent with no nuclear neon enhancement. 
For the single WN star for which neon and sulfur were 
both observed, the Ne/S ratio, which is insensitive to un- 
certainties in the star's bulk parameters, was found to be 
consistent with cosmic values, despite larger uncertain- 
ties in the total sulfur abundance. 
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